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Abstract. High-field magnetic white dwarfs have been long suspected to be the 
result of stellar mergers. However, the nature of the coalescing stars and the precise 
mechanism that produces the magnetic field are still unknown. Here we show that the 
hot, convective, differentially rotating corona present in the outer layers of the remnant 
of the merger of two degenerate cores is able to produce magnetic fields of the required 
strength that do not decay for long timescales. We also show, using an state-of-the-art 
Monte Carlo simulator, that the expected number of high-field magnetic white dwarfs 
produced in this way is consistent with that found in the solar neighborhood. 



1. Introduction 



High-field magnetic white dwarfs have field strengths larger than 10 6 G and up to 10 9 G 



(ISchmidt et al.l 120031) . The galactic population of these white dwarfs exhibits two re- 
markable pecul iarities. The first on e is that very few of them belong to a non-interacting 
binary system dKawka et al. 2007 ), wh ile the second one is that they are more mas- 



sive than average dSilvestri et al. 20071) . Generally speaking there are two competing 



scenarios to explain the formation of these stars. One possibility is that these white 
dwarfs descend fro m single stars, so the magnetic field is a fossil of previous evolution 



( Angel et"ailll98lh . However, this scenario cannot explain why they are preferentially 



massive, and wh y they are not found in non-interacting binary systems. Recently, it has 



been suggested ( Tout et ai1 l2008; Nordh aus et al.ll20lTh that strong magnetic fields can 



be produced during a common envelope episode in a close binary system in which one 
of the components is degenerate. During this phase, spiral-in of the secondary in the 
extended envelope results in differential rotation in the convective envelope, which in 
turn produces a stellar dynamo which could give raise to strong magnetic fields. How- 
ever, it has been shown that the magnetic field produced in this way does not penetrate 
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in the white dwarf, but instead decays rapidly when the common envelope is ejected 
dPotter & Toutll2010h . 

The coalescence of double degenerate objects has been the subject of numerous 
theoretical studies during the last years. Large computational efforts have been devoted 
to study this intrinsically three-dimensional problem because it is thought to be at the 
heart of sever al interesting ast r ophysical phenomena, among which we mention Typ e 
la supernova (Webbink] 1 19841 : llben & Tutukovl Il984h . magnetars (Kin g etal.1 l200lh 
and hydroge n-deficient carbon stars and of R Corona Borealis stars (Longlan d et al. 
201 lL |2012|) . Also, the large metal abundances found around some hydrogen-rich 



white dwarfs with dusty disks around them c ould also be explained by the merger of 
a carbon-oxygen and a helium white dwarf ( Garcia-Berro et all l2007h . Finally, dur- 
ing the phase previous to the coalescence of a double white-dwarf close binary system 
gravit ational waves are emitted, and it has been shown that LISA will be able to detect 
them ( Loren-Aguilar et aT1 l2005). Here we quantitatevely show that the merger of two 
degenerate cores can also explain the presence of very high magnetic fields in som e 
white dwarfs — a scenario initially suggested by Wickramasinghe & Ferrario (I2000h . 



2. An aio dynamo in the merger remnant 



Full three-dimensional simulatio ns of the merger of two white dwarfs (Guerrero et al. 



20041 : lLoren-Aguilar et al.l l2009) have shown that the final remnant of the coalescence 



is a central white dwarf containing all the mass of the undisrupted primary. On top of 
this white dwarf there is a hot corona made of about half of the mass of the disrupted 
secondary. Finally, surrounding this structure a rapidly rotating Keplerian disk can 
be found. This disk contains nearly all the mass of the secondary which has not been 
incorporated to the hot corona, since little mass (~ 10~ 3 Mq) is ejected from the system 
during the merger process. The structure of the system is shown in Fig. 1. 

The temperature gradient in the hot corona is high, and thus it is convective. Using 
the Schwarzschild criterion we found that the inner and outer edges of the convective 
region are located at R ~ 0.012/?©, and R ~ 0.026 Rq, respectively, and that the total 
mass inside this region is ~ 0.24 Mq. Assuming energy equipartition, the resulting aa> 
dynamo produces a magnetic field B 2 /Sn « p(a>R) 2 /2. For the typical values found 
in our simulations the magnetic field is B ~ 3.2 x 10 10 G. For this mechanism to be 
efficient the dynamo must work for several convective turnovers before the energy of the 
hot corona is radiated away. The temperature of the corona is very high, and hence it is 
preferentially cooled by neutrinos. Typical luminosities are L v ~ 4.0x 10 2 Lq, while the 
total thermal energy of the corona is U ~ 8.8 x 10 48 erg. Consequently, the convective 
shell lasts for Thot ~ 1.8xl0 5 yr. The convective turnover timescale is t C01W « Hp/v COIW , 
where Hp « 2.7 x 10 8 cm is the pressure scale height and v conv « 8.0 x 10 7 cm/s is the 
convective velocity. Thus, r conv ~ 3.3 s, and during the lifetime of the hot corona the 
number of convective cycles is sufficiently large. Thus, the aa> mechanism is able to 
produce strongs magnetic fields. 

Now that we have proved that large magnetic fields can be established in the af- 
termath of a double degenerate merger, it is necessary to assess if during the lifetime of 
the resulting white dwarf the magnetic field diffuses outwards, to the surrounding disk, 
or inwards, to the degenerate primary. Thus, we solved the diffusion equation for both 
the disk and the degenerate core, and we found that the timescale for diffusion of the 
magnetic field across the disk is ~ 2.0 x 10 11 yr, while for the diffusion timescale 
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Figure 1 . Temperature (solid line, left scale) and rotational velocity (dashed line, 
right scale) stratifications in the final remnant of a double white dwarf merger in 
which the binary system is composed of two stars of 0.6 and 0.8 Mq. The central 
spinning white dwarf rotates as a rigid body with a rotational velocity u> ~ 0.26 s _1 , 
and the corona rotates differentially, with a peak rotation rate co ~ 0.33 s _1 . These 
velocities arise from energy and angular momentum conservation, since little mass 
is ejected from the system, so the orbital angular momentum of the binary system is 
invested in spinning up the remnant, while the available energy is primarily invested 
in heating the corona. The location of the convective region is displayed by the 
shaded area. 



for the central white dwarf turns out to be twd ~ 4.3 x 10 9 yr. Thus, it can safely 
stated that as the white dwarf cools, the magnetic field is not allowed to diffuse across 
the keplerian disk nor to penetrate in the degenerate core, and consequently remains 
confined to the surface layers. 

Our model predicts that the masses of high-field magnetic white dwarfs should be 
larger than the average of field white dwarfs and that they should be observed as sin- 
gle white dwarfs, as observationally found . However, high-field magnetic wh ite dwarfs 
are generally found to be slow rotators dWickra masingh e~& Ferrariol |2000|) . Never- 
theless, this apparent shortcoming can be easily solved. If the rotation and magnetic 
axes are misaligned, magneto-dipole radi ation rapidly spins dow n the white dwarf. 
The evolution of the rotational velocity dBenacquista et al.l 120031) is given by to - 
-2ju 2 o> 3 sin a/(3/c 3 ), where / is the moment of inertia of the white dwarf, a is the 
angle between the magnetic and rotation axes and /u = BR^L D . Adopting typical val- 
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ues resulting from our SPH simulations we obtain a spin-down timescale t mdr ~ 
2.4 x 10 8 / sin 2 a yr, when a field strength B = 10 7 G is adopted. Hence, if both axes 
are perfectly aligned the remnant of the coalescence will be a high-field, rapidly ro- 
tating magnetic white dwarf. On the contrary, if both axes are nearly perpendicular 
magneto-dipole radiation efficiently brakes the remnant. Consequently, very young, 
hot, ultramassive, slowly rotating high-field magnetic white dwarfs can also be easily 
accommodated in our model. 

On the other hand, spectro-polarimetric observations show that in most cases high- 
field magnetic white dwarfs have fields with both toroidal and poloidal components. 
Our scenario can also account for this observational fact, since in the aco mechanism, 
convection is responsible for the generation of poloidal fields, whereas rotation is re- 
sponsible for the generation of toroidal fields. In particular, the energy available to 
generate the poloidal field component is pv 2 onv /2 ~ 4.0 x 10 20 erg, which is ~ 10% 
of the energy available to build the toroidal component, p(u>R) 2 /2 ~ 5.5 x 10 21 erg. 
Thus, we expect that the magnetic field geometry of the remnant of the merger will be 
complex. 



3. The number of high-field magnetic white dwarfs in the solar neighborhood 

To assess the number of mergers that c ould occur in the solar neighborhood we ex- 
panded an existing Mo nte Carlo code dGarcia-Berro et al" 1 119991 : fibres et all 120021 : 
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Garcia-Berro et all 120041) designed to study the galactic population of single white 
dwarfs to deal with th at of double degenerates. The population synthesis code is de- 
scribed in full depth in iGarcia-Berro et all (120121) . where all the necessary ingredients 
are detailed. Here, for the sake of conciseness, we only present the most significant 
results of our calculations. 

Our population synthesis calculations predict that the fraction of merged double 
degenerate cores in the solar neighborhood is ~ 2.9% of the total population. This 
number includes not only white dwarf mergers (~ 0.3%), but also the coalescence 
of a white dwarf and a giant star with a degenerate core (~ 1.1%), and the merger 
of two giants with degenerate cores (~ 1.5%). In these two last cases the mergers 
obviously occur during the common envelope phase, while for the cases in which two 
white dwarfs coalesce gravitational wave radiation is the final driver of the merger. 
Finally, it is important to mention as well that the distribution of remnant masses is 
nearly flat, in agreement with th e observed distribution of masses of magnetic white 
dwarfs (Nalezyty & Madejl l2004l) . 



Within 20 pc of th e Sun there are 122 white dwarfs (Holbe rget al.ll20 08). and con 



tains 14 magnetic stars dKawka et alJl2007h . Although scarce, this sample is 80% com 



plete, and allows reliably determine the true incidence of magnetism in white dwarfs. 
Mass determinations are available for 121 of them. Of the 14 magnetic white dwarfs in 
local sample, 8 have magnetic fields larger than 10 7 G, and 3 have masses larger than 
0.8 Mq — a value which is ~ 2.5cr away from the average mass of field white dwarfs. 
This mass cut is chosen to cull only white dwarfs that are expected to be the result of 
stellar mergers. Actually, our calculations predict that ~ 4 white dwarfs are the result of 
double degenerate mergers, and have masses larger than 0.8 Mq, in agreement with ob- 
servations. Additionally, our simulations predict that the fraction of white dwarfs more 
massive than ~ 0.8 Mq resulting from single stellar evolution is ~ 10%. Consequently, 
the expected number of massive white dwarfs in the local sample should be ~ 12. In- 
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stead, the local sample contains 20, pointing towards a considerable excess of massive 
white dwarfs, which could be the progeny of mergers. The rest of the population of 
magnetic white dwarfs (~ 5) — those with s mall magnetic field strengths — would be 
well the result of the evolution of single stars ( Aznar Cuadrado et al.l l2004). 



4. Conclusions 

It has been shown that the hot, convective, differentially rotating corona predicted by 
detailed Smoothed Particle Hydrodynamics simulations of the coalescence of two de- 
generate stellar cores can produce very high magnetic fields. We have also shown that 
these magnetic fields are confined to the outer layers of the remnant of the coalescence, 
and do not propagate neither to the interior of the white dwarf or to the debris region. 
Our model has two clear advantages, which meet the requirements imposed by obser- 
vations. The first one is that it naturally predicts that high-field magnetic white dwarfs 
should preferentially have high masses. The second advantage is that it explains why 
high-field magnetic white dwarfs are found to be single stars. Moreover, our scenario 
predicts that if the rotation and magnetic axes are not aligned the magnetic white dwarf 
rapidly spins down due to the emission of magnetodipole radiation. Thus, we expect 
that high-field magnetic white dwarfs should have different rotation periods, depending 
on their evolutionary status and on the geometry of the magnetic field. Moreover, in 
the case in which the masses of the merging objects white dwarfs are not equal the ke- 
plerian disk can eventually form a second-generation planetary system. Disruption of 
small bodies in this planetary system could contaminate the atmospheric layers of the 
magnetic white dwarf. This finding is important, as it can explain the recently discov- 
ered population of metallic magnetic white dwarfs. If, on the contrary, the masses of the 
merging white dwarfs are similar the remnant has spherical symmetry and rotates very 
rapidly, as observed in some high-field magnetic white dwarfs. Also, the geometry of 
the surface magnetic fields can be well explained by our model. Finally, we have also 
shown that the expected number of double degenerate mergers is roughly consistent 
with the number of high-field magnetic white dwarfs in the local sample. In summary, 
our calculations indicate that a sizable fraction of all high-field magnetic white dwarfs 
co uld be the result of double d egenerate mergers, a hypothesis previously anticipated 
by lWickram asingh e &Ferrariol ( 20001) . but not hitherto demonstrated by a quantitative 



analysis. 
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